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Abstract 



/\ • After a very brief review of the formalism of lattice gauge theories we show 

j^ ■ how one can calculate the parameters of the continuum chiral Lagrangians 

proceeding through the derivation of an effective lattice chiral Lagrangian 
as an intermediary step. The derivation is done in the strong coupling 
limit. We also discuss how the derivation could be carried out in the 
intermediate coupling domain by numerical simulation techniques. 

1 Introduction 

The lattice regularization of a quantum gauge field theory^, coupled with numer- 
ical simulation methods, provides a very powerful tool to calculate non-perturbative 
observables otherwise inaccessible to more conventional analytic techniques. Since the 
introduction of lattice gauge theory (LGT) simulations in 1979^, enormous progress 



has been made in the apphcation of these techniques. The proceedings of the inter- 
national symposia on lattice gauge theories, which have been held yearly since 1983, 
offer one of the best references for recent advances in the field^. In the course of 
the years, LGT calculations have increased in scope, scale and degree of sophistica- 
tion, progressing from the early estimates of the simplest observables to reasonably 
accurate determinations of several quantities of phenomenological interest. One of 
the features which characterizes the maturity of LGT simulations is that they have 
gone way beyond being a (poor) substitute for the laboratory provided nature, good 
only for a direct approximation of the observables of the continuum theory. Rather, 
the full power of the numerical laboratory is being more and more exploited. Thus, 
information about the dynamics of the interactions might be obtained, for instance, 
from the way in which the finite size of the lattice affects the results of a simula- 
tion. While the ultimate goal is, of course, always to use the lattice regularization 
for deriving the properties of the continuum, an increased interest in the dynamics of 
the discretized system pays off dividends with the information it indirectly provides 
about its continuum counterpart. 

In this lecture we would like to illustrate how one may derive the parameters 
of the effective chiral Lagrangians^'^'^ used for a phenomenological description of low 
energy QCD interactions from first principles by lattice techniques. The strategy 
will be to exploit the fact that the lattice regularized system also exhibits a chiral 
symmetry in order to perform, directly on the lattice, the same reduction of degrees 
of freedom which conceptually characterizes the derivation of an effective Lagrangian 
in the continuum. By taking the limit of low lattice momenta in the lattice chiral 
Lagrangian one can then calculate the parameters of the continuum chiral Lagrangian 
itself. The derivation of the lattice chiral Lagrangian will be done here in the strong 
coupling approximation. The material which we are presenting is taken from research 
which is currently in progress and a corresponding determination of the lattice chiral 
Lagrangian for more realistic values of the lattice coupling constant, which must be 
done numerically, has not been carried out yet. The strong coupling calculation, 
however, illustrates even better than a numerical simulation the reduction of degrees 
of freedom that leads to the lattice chiral Lagrangian and which is at the root of the 
method we propose. A discussion of how the same reduction could be implemented 
numerically in the intermediate coupling domain will be presented at the end of this 
paper. 



l.L Lattice QCD 

We begin by reviewing in an extremely concise manner the formalism of lattice 
QCD. By defining a quantum field theory on a lattice one achieves a regularization of 
the ultraviolet divergences which is non-perturbative and gauge invariant^. In general, 
Euclidean space-time is discretized by the introduction of a regular lattice of points 



and oriented links joining neighboring points. In the vast majority of apphcations this 
lattice is a hypercubical lattice, however, for reasons that will become clear later, we 
will need to work with a different lattice to carry out the strong coupling expansion. 
Therefore we will leave the detailed geometry of the lattice unspecified for now, but 
we will assume that the lattice is regular and uniform, with lattice spacing a. 

We will denote by x the coordinates of the lattice points and by v the fundamen- 
tal displacement vectors in the lattice. Thus, for clarification, with a hypercubical 
lattice we would have x = {iia,i2a,isa,i4^a) with integer ii,i2,i3,ii and v would 
range over (a, 0, 0, 0), (0, a, 0, 0), (0, 0, a, 0), (0, 0, 0, a). The fundamental variables are 
the gluonQ and the quark fields. 

The gauge part of the action is constructed from the plaquette variables, i.e. 
the transport factors around the elementary closed contours of the lattice. The 
plaquettes are defined by two displacement vectors, and so we will use the nota- 
tion Gx,viV2 to denote a plaquette variable. (For the hypercubical lattice Gx,vn;2 = 
Gx,viGx+vi,v2Gx-^v2,viGl: ^2)- ^^^ gS'Uge part of the action^ is then given by: 

S, = PK, J2 Tr(/-G,,,,,J, (1) 

the coupling parameter /? being related to the bare coupling constant by /? = G/g^. 
The plaquette variables generalize the field strength tensor F^^{x) of the continuum 
theory and the factors in Eq. 1 are arranged so that the lattice action reduces formally 
to the standard action of the gauge field in the limit of vanishing lattice spacing 
a — > if the gauge variables are identified with continuum transport factors through 
Gx,v = exp[igoA^{x)v^']. 

The quark fields are defined over the sites of the lattice and are given by elements 
of a Grassmann algebra ipx (color, spin and flavor indices will frequently be left 
implicit). The quark matter-field action is given by: 

Sg = K2J2 J2('^^'^^^,^^^+'" ~ i^x+vi^Gl^^ipx) + Yl i'xjxi^x- (2) 

XV x 

In terms of fundamental variables the (lattice regularized) quantum expectation 
value of any observable 0{G, ijj, tp) is given by 

<0>= Z-' J DGx,. DijxD^xO{G, ^, ij) exp [-Sg{G) - S,{G, ^, ^)] (3) 



with 



Z = JdGx,v D^xD^exp [-S,iG) - Sg{G,^,^)]. (4) 



'^ Gauge fields are customarily denoted by Ux^v, but eventually, following another standard nota- 
tion, we will want to reserve the U symbol for the chiral fields. We will therefore denote the lattice 
gauge fields by Gx.v 

^Ki and K2 are suitable normalization constants, dependent on the actual shape of the lattice. 



It is to be noticed that, if one considers first a lattice of finite volume V, letting 
l^ ^ oo only at the end of the calculations, all of the integrals in Eqs. 3 and 4 are 
finite and gauge invariant. Thus the lattice formalism provides a mathematically well 
defined, non-perturbative and gauge invariant regularization of QCD. Formally, in 
the limit of vanishing lattice spacing, the above expression for the quantum averages 
of the observables reduces to the continuum path integral, but, of course, the latter is 
per se not well defined. The correct continuum limit is obtained through the process 
of renormalization, whereby lattice spacing a and bare coupling constant go are sent 
simultaneously to according to a precise functional relationship which, because of 
the property of QCD of being an asymptotically free theory, can be established by 
perturbative techniques. 



1.2. Effective chiral Lagrangians 

The most important property of QCD for low energy phenomenology is that the 
theory has an approximate S'f/(3) x S'f/(3) chiral symmetry, which would be exact 
for niu = iTLd = nis = 0, and that this symmetry is spontaneously broken. As a 
consequence, the spectrum of particle states contains eight pseudoscalar mesons of 
very low mass, also called pseudo Goldstone bosons, which would be the Goldstone 
bosons of the spontaneously broken symmetry if the masses of the lightest quarks 
were indeed equal to zero. 

Effective chiral Lagrangians have been considered in other lectures delivered at 
the ICTP summer workshops and we refer the reader to the proceedings or to the 
seminal papers of Gasser and Leutwyler for more detailed information. We summarize 
here only the most basic features of effective chiral Lagrangians and of their use for 
low energy QCD. 

Following Gasser and Leutwyler^'^, one introduces the generating functional for 
connected Green's functions of QCD in the presence of scalar, pseudoscalar, vector 
and axial vector sources 

Z{s,p,v,a) = /'DgDgDG'e'/'^'"^^^'^"^'^'^'"''^), (5) 

^ = ^QCD - Qi^) [s(a;) - i-f5P{x)]q{x) + q{x)Y[v^{x) + 75a^(x)]g(x). (6) 

In this expression, q, q represent the lightest quark flavor triplet and CPqcd i^ 
what remains of the Lagrangian of QCD when the masses of three lightest quarks are 
set to zero 

^QCD = -^,G,.{x)G^^{x) + q{xh^[td, + G,{xMx). (7) 

s{x),p{x),Vfj,{x) and a^lx) are Hermitian matrices in color space. The effective chiral 
Lagrangian is introduced to provide a description of phenomena involving small ex- 
ternal momenta. Such phenomena will lead to the excitation of the pseudo Goldstone 



bosons only, and correspondingly the cliiral Lagrangian will be formulated in terms 
of an octet of pseudoscalar fields 7r'^(x), collected in a unitary 3x3 matrix: 

Uix)=exp[^^^^^^], Uix)Uix)^ = l. (8) 

A", a = 1, ..., 8 are Hermitian generators of SU{3), with the normalization 

Tr(A''A^) = 25"'', (9) 

F^ is the pion decay constant in the chiral limit: 

A" 

< I g(a;)7M75Yg(a;) | tt^ > = iUPf,6''\ 

U = Fo{l + 0{mg)) ^ 93.3MeV. (10) 

The chiral Lagrangian itself takes the form of a perturbative expansion^ in the 
external momenta p: 

Z{s,p,V,a)= /■/}f/e^/'^'^^2(C/;s,p,«,a)+£4(C/;s,P,^,a)+...^ ^^^) 

Here £2 is a Lagrangian involving terms of order p^: 

£2 = ^Ti{[D,Uix)][D^^Uix)]^ + x{x)Uixy + xix^Uix)}. (12) 

In this expression, 

D^U{x) = d^U{x) - ir^{x)U{x) + iU{x)l^{x) (13) 

is a flavor covariant derivative and matrix x{^) collects scalar and pseudoscalar 
sources: 

X{x) = 2Bq[s{x) + ip{x)]. (14) 

The constant Bq is related to the quark condensate in the chiral limit: 

< I MM I >o=< I Jrf I >o=< I ss I >o= -F^Bo{l + 0{mq)). (15) 

Similarly, £4 is the most general Lagrangian of order p^: 

£4 = L,{i:i[{D^U){D^U)^]Y + L^i:i[{D^Uy{D,U)\Vi[{D^U)\D^U)] 

+ L^i:i[{D^U)\D>^U){D,U)\D''U)] + L4Tr[(D^[/)t(D'^f/)]Tr(x^f/ + xU^) 

+ L,i:i[{D,U)\D^U){x^U + xU^)] + HTiix^U + xU^)f 

+ Lj[Ti{x^U - xU^)f + L^-TtU^Ux^U + xf/^xf/^) 

- iL^Ti[F^,{D^U){D''U)^ + F^^,(D'^?7)t(D"f/)] + LioTr(t/^F^^,f/F^'^'') 

+ giTr(F^^,F^^'- + F^^.F^ ^'-) + //2Tr(xtx)- (16) 

'^It can be shown that in this expansion vector and axial vector sources are of order p, whereas 
scalar and pseudoscalar sources are of order p^ . 



Here F^ , Fj^'^ are field strengths for riglit and left-lianded fields: 



Fr = d^'r'' - d'^r" 



i\r^, r 



R - L/ / - L/ / - t^/ , / J, 

As we see from the above, at leading order, chiral symmetry restricts the number 
of terms in the effective chiral Lagrangian to only two, so we need to know only two 
constants, Fq and Bq, in order to describe all low energy phenomena in QCD. At next 
to leading order, we need ten more constants, Li - Liq. (Constants Hi and H2 are of 
no physical significance.) Typically, these coupling constants are fixed by comparison 
of the predictions that follow from the Lagrangian with a subset of the experimental 
data. One can then derive further predictions that can be tested against experiment 
or used as input in the study of other interactions (e.g. one might need to calculate 
strong corrections to electroweak matrix elements). 

Of course, it would be very important to be able to calculate the values of the 
coupling constants appearing in the chiral Lagrangian from first principles, directly 
from the fundamental QCD Lagrangian. This requires that one carries out explicitly 
an integration over the high energy degrees of freedom of QCD, trading off quark and 
gluon fields for the fields of the pseudo Goldstone bosons. In the phenomenological 
applications one assumes the resulting form of the Lagrangian, but, as we have just 
stated above, one derives its parameters from comparison with experiment. What 
we will show in the remainder of this lecture is that the integration over quark and 
gluon fields can actually be performed explicitly in the lattice formulation in the 
large N, strong coupling limit and that it leads to a lattice chiral Lagrangian from 
which the coupling constants of Eq. |16| can be derived. Finally we will discuss how 
it may be possible to use numerical simulations to go beyond the strong coupling 
approximation and derive the chiral Lagrangian for more physical values of the lattice 
coupling constant go. 

2 Lattice chiral Lagrangian in the strong couphng 
approximation 

In order to derive a lattice chiral Lagrangian one must integrate out the high 
energy degrees of freedom, i.e. the quantum fluctuations of the gauge and quark fields, 
re-expressing the generating functional in terms of pseudoscalar fields conjugate to 
the external sources. Thus one effectively performs a bosonization of the original 
theory. For general values of the bare lattice coupling constant this can only be 
done numerically, but in the strong coupling and in the large N (number of colors) 
limits, the bosonization can be done analytically. The techniques for deriving an 
effective Lagrangian in the strong coupling limit were laid down in pioneering papers 
by Kluberg-Stern, Morel, Napoly and Petersson''^ and Kawamoto and Smit*^. We add, 
however, to the work of these authors an ingredient that plays a crucial role for the 
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possibility of obtaining the terms of order p^ in the chiral Lagrangian. As in most 
lattice calculations, Refs. 7 and 8 deal with a hypercubical lattice. While such a lattice 
obviously does not have the symmetry of the continuum, as is well known all tensors 
of rank up to two symmetric under the group of lattice transformations are also 
symmetric under the full group of 4-dimensional rotations of Euclidean space-time. 
However, this property does not carry through to tensors of higher order, so that, 
although one could expand the chiral Lagrangian derived on a hypercubical lattice 
to terms of 4*^^ order in lattice momentum, one would obtain terms that cannot be 
identified with continuum counterparts. 

Inspired by the use of lattices of higher symmetry in the theory of lattice gases^, 
we will formulate the theory and perform the strong coupling expansion on a body- 
centered hypercubical (BCH) lattice^. The BCH lattice can be obtained from a 
hypercubical (HC) lattice in two equivalent ways. One can define it as a HC lattice 
with all the sites of odd (or even) parity removed. Equivalently one can take as sites 
of the BCH lattice all the sites of a HC lattice together with all the centers of its 
elementary cells (hence the name BCH). In the BCH lattice every site has 24 nearest 
neighbors. As fundamental displacement vectors v in the positive direction we will 
take the following: 

< = ^^^,{l<^<J<4,a = ±l}. (18) 

The BCH lattice has the largest symmetry group of all four dimensional lattices, 
with 1152 elements (for comparison, the symmetry group of the HC lattice has 384 
elements). As a consequence of this larger symmetry group, all tensors of rank up to 4 
with the symmetry of the BCH lattice are also invariant under the full 4-dimensional 
rotation group. This property will be very important for our derivation. 

The use of the BCH lattice entails another advantage in the strong coupling 
limit. It is of rather technical nature and so we will mention it only briefly. The 
definition of lattice fermions encounters some notorious problems. A straightforward 
discretization of the Dirac operator leads to the introduction of additional poles in the 
propagators of the fermions at the corners of the Brillouin zone (species doubling). 
Although these extra fermionic modes would in general be coupled through the gluonic 
quantum fluctuations, on a hypercubical lattice there is a fourfold degeneracy of 
totally decoupled fermionic degrees of freedom. This is akin to starting with one 
flavor of fermions, but discovering that the theory actually contains four decoupled 
flavors. As a consequence, the spontaneous breaking of chiral symmetry gives origin to 
16 pseudo Goldstone bosons, even though one started with a single flavor. On a BCH 
lattice, because of the higher coordination number of its vertices, such a separation 
of the fermionic flelds into four decoupled components (which is not demanded by 

'^ This lattice, also known as F/^ lattice, has been considered in prior unrelated LGT investigations, 
see for instance Refs. 10-13. However it is the motivation for its use in the theory of lattice gases 
that comes closer to the one underlying our own work. 



any general theorem on lattice fermions) does not take place and one maintains the 
relationship between number of explicit flavor indices and number of pseudo Goldstone 
bosons proper of the continuum theory.^ 

Our goal is to express the generating functional of Eq. 5, formulated over the 
BCH lattice (cf. Eqs. 1-4), in terms of lattice bosonic variables [4 = exp [in'^X"- / F^. 
We follow the work of Refs. 7, 8, adapting it to the fact that we are on the BCH 
latticeQ. Also, for the sake of brevity, we will not include into our equations explicit 
reference to the external sources v and a coupled to the vector and axial vector 
currents. These can be introduced via suitable chiral transport factors defined over 
the links of the lattice. We will include the results for the relevant corresponding 
coefficients in £4 only in the final table. 

In the strong coupling expansion, the gauge term (Eq. |l|) is suppressed by 1/ Qq- 
Therefore, we can neglect it to the leading order in the strong coupling expansion, 
and the action reduces to Eq. ^. Now the gauge variables appear only in the first two 
terms, moreover they appear only linearly. This means that the product of integrals 
over gauge variables factors into a product of one-link integrals: 

Z= I D^jD^e- S- "^^^^^^ n / dGx,ve~^^'^ ■ (19) 

•^ x,v '' 

Here we have defined: 

W^,^ = Y^(V^x?^G'^,i,^x+t; - i^x+vfGl^^ip^). (20) 

One link integrals of this type have been calculated for the case when the gauge 
field couples to a bosonic source in Ref. 14. In Ref. 7, it was shown that this result 
applies equally to an interaction of gauge fields with fermionic fields. The interested 
reader should consult Ref. 7 for the details. Here we just give the final result for the 
generating functional to the leading order in the strong coupling, large-N limit: 

Z= f D4jD4jexp{NY,TT[F{X{x,v))]-Y,tP,j,tP,}. (21) 

X,V X 

In this equation, \{x,v) is a matrix in Dirac-fiavor space, defined by its elements^: 

1 N An 



'^It has been argued in Ref. 11 that the use of the BCH lattice actually worsens the problem 
of lattice fermions as one tries to recover the continuum limit. This is irrelevant for our purposes, 
because our motivation for its use comes solely from the need of obtaining better symmetry properties 
in the strong coupling limit. We would not advocate using the BCH lattice in numerical simulations 
done in the scaling regime, where there is good evidence for the recovery of rotational symmetry 
with the ordinary HC lattice. 

•''On this lattice the values of constants in Eqs. El and are: Ki — K2 — 1/12. 

^ Greek letters will denote indices in the direct product of Dirac and flavor spaces, and Latin 
letters the color indices. In this section we will take the number of flavors to be n and will specialize 
to n = 3 only at the end. 



F{X) is the following matrix function of A: 



l + v^r^A 



F(A) = l-Vr^A + ln^ , (23) 

which should be truncated for finite A^ as a power series in Grassniann variables. 
Therefore, Eq. ^is valid only in the large-N limit. 

Next we proceed to convert the fermionic path integral into one over bosonic 
variables as in Refs. 7, 8. In Eq. ^, A^Tr[F(A)] is the interaction term. It is important 
to notice that it depends only on color singlet combinations ipf {x)il!f (x) at different 
points on the lattice. Since these appear also in the source term, —Y^xi^xjxi^x, we 
can rewrite the full generating functional as: 

Z = e^S...^[^(^(^'^))l y'D^Z}V^e-S.'^"^>'^% (24) 

with: 

1 ^" 9 d 

A(...)». = ^, EJgJ^ifh^g^yU (25) 

Now we have to perform the integration over fermionic variables. This integral splits 
into a product of one-site integrals: 

Zo= f DijDije-'^ = n / dij^d^^e-^-^-"^- = \{ zo{jx), (26) 

•^ X "^ X 

zo{jx)= fdiJxd^xC-^^'^^^. (27) 



There are two methods of integration over the fermion fields. The first method^ is 
based on the use of Laplace transforms. We will follow the method of Ref. 8, where 
it was shown that ^o(jx) can be written in the following form: 

zoUx) = /'rfA^^e^^(^'^^-)-^^(''^^^)+™"^*. (28) 

Here A^^ is a unitary bosonic matrix in Dirac-flavor space, the trace is over Dirac- 
flavor indices and the constant is irrelevant so it will be omitted hereafter. 

Using Eqs. ^ and ^ we obtain an expression for the full generating functional: 

Z= f DMe^^^^.^"^^^^^''-^^^^^-"^^'''^''^-^^^'^^''^''^^ = j DMe^^^\ (29) 

where A is now a bosonic matrix: 

X{x,v) = UMx+.fMx. (30) 

y 



We look for a translationally invariant saddle point of the action in (ESf) of the form 



Ml, = V, (31) 
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V being proportional to the unit matrix in Dirac-flavor space. From this we obtain 
the effective potential: 

K =J-^^^. (32) 

volume 



The saddle point for the massless case is given by: v{am = 0) = V23/4 ^ 1.20. We 
will parameterize the field TVIj. in a nonlinear way*^: 

M. = vexp [^{S, + 7r,.75 + ^^7. + ^^^7.75 + ^T^a,,)]. (33) 

Here S,7r,V'^,A'^ and T'^'' are sixteen Hermitian matrices in the flavor space. By 
finding the propagators of these fields, one can show that only vr are Goldstone bosons. 
In principle, we would need to integrate out the scalar, vector, axial vector 
and tensor fields in order to obtain an effective action for pseudoscalar mesons. In 
first approximation we will neglect these contributions^, and concentrate only on the 
direct interactions among Goldstone boson^: 

M. = vexp{z^). (34) 

-To 

We will expand the action in Eq. ^ up to the fourth order in Taylor series 
around the vacuum. This will be sufficient to extract coefficients of effective chiral 
Lagrangian to 0{p'^) since the deviation from the vacuum A^,^ — Aq is of 0{p). 

The action in Eq. ^ can be written as: 



^ = ^EE{rTUvr Tr[(A.,.-Ao)1} + iV^ETr(j.A^. 

x,v n=l "'• \ "^^ / Ao ^ 



(35) 



We will evaluate the Dirac trace in a basis where 75 = diag(l, 1,-1,-1). By 
using the relation 

^M.+, = Ml^J, (36) 

which follows from the definition oi Ai, we can take the Dirac traces to obtain: 

Tr[(A,,, - Ao)"] = 2X^TTf[{uU,U, - 1)" + h.c. ]. (37) 

Here we have defined new variables which are matrices in fiavor space only: 

[/, = exp(z^). (38) 

-TO 

''They are suppressed by the masses of relevant resonances. 

*In this paper we will not consider the effects of the U(l) anomaly in any detail, so we will restrict 
our attention to the traceless part of tt^ ■ 
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Try indicates a trace over flavor indices, and from now on we will omit the subscript 
/. Evaluating the second term in Eq. ^ is straightforward and one obtains the 
expression for the action in terms of new variables: 

S = SiVEE-rh^ Tr[(t/U[/.-ir+h.c.] 

+ 2iVt;ETr(?«t^x+ h.c. )• (39) 

X 

From now on we will denote by a subscript v any variable which under the lattice 
symmetry transforms in the same way as the link f °- on which it is defined. For 
example: 

a^ = 1^-^, a^ = (a^,) , etc. (40) 

V2 

Partial derivatives in the direction of the link f ^^ are defined similarly: 



«j 



ditt + adjtt ^2 



dya = _ , d^a = dy{dya), etc. (41) 

v2 

These derivative^ appear in the Taylor series expansion of Ux+v in powers of a: 

Ux+y = U + a{dyU) + ^(d^yU) + . . . , (42) 

One can now substitute the previous Taylor series expansions into Eq. |39|. Using the 
following results for summation of tensors of rank two and four (which are easy to 
prove from definitions in Eq. ^0|) over the twelve links in the positive direction: 

EaA = 3E«A«&^ (43) 

V fl 

E CLvKcvdy = - Y^{a^h^'aulf + a^a^lf + a^Ka^'h^), (44) 

one obtains a manifestly Lorentz invariant expressions for the action. We will not 
give the details^ of this rather straightforward but tedious calculation and just give 
the final result in the following table. 

^Hcre we are using the convention that, when there is no subscript x, a variable (or its derivative), 
is evaluated at x. 

*^The interested reader should consult our longer paper, currently in preparation, where we give 
all the details as well as a discussion on how one includes currents. 
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Table 1: Comparison of the strong coupling expansion with experimental results from 
Refs. 6, 18. The left column contains our results from the strong coupling expansion. 



The right column contains the experimental values L[ 
of 10^3^ 



.rrir^ 



The numbers are in units 





theory 


experiment 


Li 


0.5 {Foaf 


0.9 


± 0.3 


L2 


1.1 {Foay 


1.7 


± 0.7 


L, 


10.4 (Foa)2 


-4.4 


± 2.5 


U 


(Foaf 





± 0.5 


U 


(Foaf 


2.2 


± 0.5 


U 


(Foaf 





± 0.3 


Lr 


0.9 (Foa)2 


-0.4 


± 0.15 


Ls 


2.6 {Foaf 


1.1 


± 0.3 


L9 


-6.7 {Foaf 


7.4 


± 0.2 


Lio 


-6.7 (Foa)2 


-5.7 


± 0.3 



3 Discussion 

We have chosen to express the entries in Table 1. in terms of a common factor 
(Foo)^ because of a point we will be making below, but the value of this quantity is 
also determined by the strong coupling expansion, from the O(p^) terms in the lattice 
chiral Lagrangian, and is given by 2.1A^. With 3 colors the entries in the table come 
out about one order of magnitude higher than the values derived from experiment. 
We take this to be an indication that the strong coupling limit produces too tight a 
binding between quark and antiquark, which in turn leads to an unacceptably high 
value when Fq is expressed in terms of a~^ (or, equivalently, to too large a value for 
a). If we allowed ourselves to set a^^ ~ Fq, which would be a more reasonable scale 
for a strong coupling calculation, then the theoretical predictions for several of the 
coefficients would compare rather well with the experimental results. Of course we 
should not make too much of this agreement, because we expect the values of the 
coefficients to depend on detailed dynamical features of the interactions which may 
not be well reproduced by the strong coupling approximation. What is much more 
important is the demonstration, through the strong coupling expansion, of the main 
point we wanted to make, namely that one can formulate a chiral effective Lagrangian 
on the lattice and that this can be the vehicle for the derivation of the parameters of 
the continuum chiral Lagrangian. 

Another point that should be mentioned is that the expansion for small lattice 
momenta which we have used to derive the coefficients corresponds to using the 
lattice chiral Lagrangian in the tree approximation. We have calculated a few of the 
corrections that would be induced by one loop diagrams and have found these to be 
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small. However, the same factor of (-Fo^)^ which we discussed above, enters also as a 
coefficient in the denominator of the loop diagrams, so that the statement that the 
loop corrections are small should be taken with caution. They would become larger 
if one could assign to F^a the smaller value which we advocated above. However, 
even if the loop corrections turned out to be large in a more realistic calculation, 
this would only represent a technical problem and not a conceptual difficulty for 
the whole approach. The very important fact is that, because we are working with 
a regularized theory, all physical quantities are well defined and finite. Moreover, 
since the lattice system also exhibits spontaneous breaking of chiral symmetry and 
Goldstone bosons in the limit of vanishing quark mass, it will always be possible to 
perform an expansion for small lattice momenta. The coefficients of such expansion 
are in any case well defined quantities. Thus the question is whether they can be 
calculated by a perturbative expansion of the chiral effective theory, which would 
be more convenient, or whether, failing such possibility, one will have to resort to 
numerical techniques. Nevertheless, even in the latter case, one would still be dealing 
with a bosonic system, and would therefore avoid the need of simulating the quantum 
fluctuations of fermionic variables, which represents today the major difficulty one 
faces in the implementation of QCD numerical simulations. 

Finally, in order to obtain reliable values for the coefficients one should perform 
the calculation in the range of values of the lattice coupling constant {j3 = Q/qq ~ 
5.7 — 6, including quark degrees of freedom in the simulation) where one witnesses 
the onset of scaling to the continuum. Such a calculation can at present only be done 
by numerical techniques. The way we envisage it could be carried out would be by 
assuming a sufficiently large set of couplings for the lattice effective theory and fixing 
them through the matching of an overcomplete set of expectation values. This is very 
similar to procedures commonly used in numerical studies of renormalization group 
transformations. The crucial point is that, because the lattice effective theory already 
accounts for the long range excitations of QCD, the matching should only require a 
reasonably small lattice size, of the order of the inverse p mass. The two theories 
(the one derived from the original QCD Lagrangian and the effective theory) should 
produce exactly the same values for the observables on any lattice size, because they 
are mathematically equivalent. Working with a small lattice, we are confident, one 
would be able to fix the parameters of the lattice effective theory with a good degree 
of accuracy and would thus establish a solid base for the subsequent determination 
of the parameters of the continuum chiral effective Lagrangian. 
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